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ASYMPTOTIC DETERMINANT OF DISCRETE 
LAPLACE-BELTRAMI OPERATORS 

ANANTH SRIDHAR 


Abstract. We study combinatorial Laplacians on rectangular subgraphs of d? 
that approximate Laplace-Beltrami operators of Riemannian metrics as e —>■ 0. 
These laplacians arise as follows: we define the notion of a Riemmanian metric 
structure on a graph. We then define combinatorial free field theories and describe 
how these can be regarded as finite dimensional approximations of scalar field 
theory. We focus on the Gaussian field theory on rectangular subgraphs of 1? 
and study its partition function by computing the asymptotic determinant of the 
discrete laplacian. 


1. Introduction 

On a graph X = (V, E), the graph laplacian is an operator that acts on the space 
of fnnctions / : E —)■ M as: 

N{u)= ^ f{v)-f{u) (1.1) 

[u,v]^E 

The determinants of graph laplacians have been well stndied in many different con¬ 
texts. For instance, in spectral graph theory, Kirchoff’s theorem relates the de¬ 
terminant of the graph laplacian on X to the nnmber of rooted spanning trees on 
X. Similarly, by Temper ley’s trick and its generalizations, the determinant of the 
graph laplacian on certain graphs X compntes the number of perfect matchings on 
a related graph X' [KPWOO]. In each case, the determinant of the laplacian is the 
partition function of an associated model in statistical mechanics. 

When X is the graph of a regular lattice, the asymptotics of the determinant have 
been investigated in many papers. In particular, when A is a rectangular subgraph 
of with NxM vertices, the determinant of the laplacian (with Dirichlet boundary 
conditions) as A, M —)■ cxo is computed by [DD] in the context of Hamiltonian walks 
on the grid. The asymptotic expansion is: 

AC 1 

logdet Xmm = NM — + 2(A + M) log(\/2 - 1) + - log(AM) + 0(1) (1.2) 

71 A 

where G is the Catalan constant. More generally, the asymptotics of the deter¬ 
minants of laplacians on rectilinear regions of have been computed in [KenOO]. 
When appropriately rescaled, this graph laplacian can be understood as a hnite 
difference approximation of the laplacian ^ ^ of Eulidean space. 

In this paper, we study a generalized graph laplacians in Z^ that can be regarded 
as hnite difference approximations of Laplace-Beltrami operators induced by a Rie¬ 
mannian metric g. These laplacians arise as follows: on any graph X, the cochain 
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complex (C*(X), d) can be regarded as a combinatorial version of the de Rham com¬ 
plex. We define a Riemannian structure on X as a set of inner products (, )c*{x) on 
the cochain spaces. On such a graph with Riemannian structure, we dehne a combi¬ 
natorial Guassian held theories, as outlined in section §2.3. The partition function 
of the boson held theory is Z = det A, where the combinatorial laplacian A is by 
{df,df)cHx) = (A/,/)c*(x) for all / G C^X). 

Assume now that the underlying space M C of A is endowed with a metric 
9 = Qxx dx^ + gyy dy'^. By square subdivision of A, we have a sequence {A*} of 
graphs with mesh —)■ 0. Assume that the {A*} have Riemmanian structure such 
that the combinatorial laplacian A* on A* is a good approximation of the laplacian 
induced by on M (see section 2.2.3). This can be regarded as discretization of a 
massless scalar held theory on a nontrivial background metric. We show that the 
determinant of A* has the asymptotic expansion (Theorem 3.1): 

logdetA* = 4 / F{g)dxdy + -[ B (g) ■ ds + + 0{1) (1.3) 

G Jm G JdM ^ 

where the functions F and B are given in equation (3.2). 

The discrete Guassian held theories are known to be conformally invariant in the 
scaling limit i —)■ oo, and consequently the asymptotic expansion of the form (1.3) 
is generally anticipated by the theory of hnite size scaling [GaPes]. In particular, 
the coefficient of the logarithmic term is universally related to the universal charge 
of the limiting conformal held theory and the Euler characteristic of the M. 

Some work has been done on relating the constant terms in the asymptotics of the 
determinant of discretized or hnite diherence laplacians with regularized functional 
determinants of analytic laplacians. In one-dimension, there are many positive re¬ 
sults, for example [For92], [BFK]. In two dimensions, for some particular cases such 
as hat rectangles and hat torii, with particular choice of discretized laplacians, the 
constant term is known to be equal to the zeta-regularized determinant of the ana¬ 
lytic laplacian, [DD],[GJK10]. We do not attempt to calculate the constant term in 
the asymptotic determinant (1.3) and its investigation for further work. 

The structure of the paper is as follows. In §2, we describe the notion of a 
Riemannian structure on a graph, and dehne the combinatorial codiherential and 
combinatorial laplacian. We explain how a Riemannian manifold can be approxi¬ 
mated by a sequence of graphs with Riemannian structure. Then we outline the 
dehnition of a combinatorial Gaussian free held theory on a graph and describe the 
relation to the dimer model of statistical physics. In §3, we focus on the bosonic held 
theory on {M,g), when M is a rectangular region in endowed with a diagonal 
metric g, and set up the calculation of the determinant. Sections §4 and §5 are more 
technical sections: §4 is devoted to estimating the inverse (A®)“^ as z —)■ oo] in §5 
we prove the main theorem, equation (1.3). Finally in §6, we conclude with some 
general remarks and further research directions. 
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regarding combinatorial Laplacians and dimer models, and for proofreading some 
versions of this text. 


2. Preliminaries 


Let X be a cell decomposition of a simply connected region M C with g-cells 
Xq for q = 0,1,2. The space of g-chains Cq{X) for g = 0,1,2 is the real vector 
spaces with canonical basis Xq. The dual space space C"^(X) = Cq{X)* is the space 
of functions Xq —)■ M. There is an isomorphism Cq{X) ~ C'^(X) induced by the 
canonical basis on Cq{X). For any subcomplex L G X, there is a natural projection 
p : C"^(X) —)■ C'^{L) given by restriction to L. Define C‘^{X,L) = ker(p) C C"^(X). 
The boundary of X is the 1-dimensional subcomplex dX contained in dM. 

On an oriented complex, there is a natural boundary operator dq : Cq{X) —)■ 
Cq-i{X). The coboundary dq : C''^(X) —)■ is the adjoint of dq with respect 

to the pairing 0^(X) x Cq{X) —)■ R. The coboundary operator satisfies dq+idq = 0 
and the complex {C*{X),d) can be regarded as a discrete analog of the de Rham 
complex of a smooth manifold. 

In the remainder, we assume for simplicity that X is embedded such that image 
of each 2-cell is a simple convex polygon. In this case the 1-skeleton of X then 
defines a simple planar graph embedded in M. 

2.0.1. The dual and double of a graph. We define the dual complex X as the complex 
found by truncating the standard dual, so that X is also a cell complex of M (see 
figure 2.0.1). Note that as usual, there is a bijection from Xq and to {X\dX) 2 -q, 
but in addition the truncation gives a bijection of dXq and dXi_q. 

The double D{X) of X is, roughly speaking, the union of the X and *X (see 
figure 2.0.1). Precisely, D{X) is the complex consisting of: a 0-cell Pa at the center 
of each a G X; a 1-cell [pa,PT] for each pair a G X*, r G Xj+i with a adjacent to 
r; and a 2-cell [pa,Pi 3 i,P'y,Pi 32 ]j for each set of a G Xq, fdi, (^2 G Xi, 7 G X 2 pairwise 
adjacent. Note that the graph of D{X) has a natural bipartite structure given by 
painting p^ white if a is a 1 -cell, and black otherwise. 



Figure 1. A graph, its dual, and its double. 


2.1. Riemannian Structures on Graphs. Recall that the metric on a compact, 
oriented Riemannian manifold (M, g) induces a scalar product on the de Rham 
complex for f,g E 



where * : f\^TpM —)■ ^TpM is the Hodge dual. Continuing to view the cochain 

spaces C^{X) as the combinatorial analog of we define: 
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Definition 2.1. A Riemannian structure on X is a set of scalar products {,)ci{x) 
for g = 0,1, 2. 

The scalar products can be required to satisfy certain additional axioms. In 
particular, recall that the scalar product on is local in the sense that it is the 

integral of a fiberwise scalar product. We consider more generally the combinatorial 
notion of a fiberwise map —)■ Q^{M) as follows: 

Definition 2.2. A map T : C'^(X) —)■ C^(X) is said to m-local if the matrix 
elements in the canonical basis of C*{X) vanishes whenever the graph theoretic 
distance between Po-,Pt ^ D{X) is greater than 2m. 

Definition 2.3. A Riemannian structne {,)c<i{x) is local if the corresponding dual 
operator C^lX) C^lX) is m-local for some m. A Riemannian structnre is said 
to be diagonal if it is 0-local. 

A diagonal Riemannian strncture can be understood as an assignment of a volume 
weight |(t| = (n, (t)c'i(x) to each cell a G Xg. 


2.1.1. Approximation theory. We now consider how a Riemannian manifold can be 
approximated by a graph with Riemannian structure. We assume as before that 
X is a cell decomposition of a region M C but in addition assume that M 
is endowed with a smooth Riemannian metric g. Define the mesh of X is ex = 
snp{ diam(r) | r G X}. The embedding of X in M defines a surjection, namely the 
de Rham integration map Ix ■ —)■ C'^(X): 


Ix(^ = “ / 

a<=X„ 


u 


for u G 

Let {X*}jgis} be a sequence of complexes in M, such that ex^ 0, for example 
generated by iterated snbdivision of an initial complex. A sequence of Riemannian 
structures on the {X*} is said to be consistent if the integration maps lx* approaches 
an isometry of spaces VL^{M) and C"^(X*) as i —)■ oo. More precisely: 


Definition 2.4. A seqnence of Riemannian structures on {X*} is consistent if for 
any / G there exists a constant pf snch that, for all i: 


{^xif 1 -^x*/)c9(Aq 


(/, f)n<i(M) 


< Vf ex* 


2.2. The Combinatorial Laplacian. Let X be a complex with a Riemannian 
structnre. Let L C dX be a subcomplex of dX. Let dg^i '■ C‘^iX,L) —)■ C^{X,L) = 
Tiid where til ■ C"^(X) —)■ C'^{X,L) is orthogonal projection. The combinatorial 
codifferential combinatorial codifferential dg^L '■ C^(X) —)■ C*^~^{X) is defined the 
adjoint of dg. 

{dq,Lf, 9 )c‘i+^{X,L) = (/, <^g+l,Lfi')c'J(X,L) 
for all f e Cl{X),ge C"'+i(X). 
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Remark 2.5. Choosing the canonical basis for the C*{X) and writing the scalar 
products (, )ci{K) as matrices Qq, the codifferential 5q can be written in matrix form 
as: 

^<1 ~ 9q-ldq9q 

From this, it is clear that generally 6q is m-local (for any ml) only if the scalar 
product Qq^i is diagonal. 

2.2.1. Duality. In this section, we briefly explain how Hodge duality arises in this 
combinatorial framework. Let X be a complex with a diagonal Riemannian struc¬ 
ture, and as above L C dX. Let L C dX be the complement of the image of L 
under the bijection dXq —)■ dXi^q, and be the complex found by removing from 
X the image of L under both bijections Xq —)■ {X\dX) 2 -q and dXg —)■ dXi_q. 

There is a bijection of Xq\Lq and X 2 -q\L 2 -q, which induces an isomorphism 
0 : Cq{X,L) —)■ C 2 -q{XL, L). Putting a Riemannian structure on Xl via the 
isomorphism, 0 becomes an isomorphism of chain complexes {Cq{X,L),di) and 
(C 2 _g(Xi„ L), 6i); ie. 0 o dq^L = i o 0 and 0 o dq^i = d^x ° 

2.2.2. The Combinatorial Laplacian. The combinatorial Laplacian is dehned as Ag= 
dq+i^Ldq^L + dL,q-idq,L- This Laplaciau should be understood as corresponding to 
mixed boundary conditions, with relative (Dirichlet) boundary conditions on L and 
absolute (Neumann) boundary conditions on dX\L. In particular, since we focus 
attention on 0-forms, we let: 

Al = 

For a diagonal Riemannian structure, it follows by straightforward computation 
that for p G Xo\L, the combinatorial laplacian can be written: 

Al/(p)= 5^ hi (/(,)_/(p)) (2.1) 

IpI 

p n (TT^o 
a n L=0 


2.2.3. Consistency of the Combinatorial Laplacian. Fix a submanifold L C 0M, 
and let A be the Laplace-Beltrami operator acting on C°° functions with Dirichlet 
boundary conditions on L and Neumann on dM\L. Consider a sequence of graphs 
{Xj} approximating {M,g) as in Section 2.1.1, and such that for each i there is a 
subcomplex Li C dXi that is a decomposition of L. In this context: 

Definition 2.6. The sequence of combinatorial laplacians A* on Xj is said to be 
consistent if for any / G C°°{M) there exists a constant rjf such that, for all i: 

Ai/x,/-/x,A/l<pp4_ (2.2) 


2.3. Gaussian Field Theories on Graphs. In this section, we briefly outline the 
dehnition of discrete Guassian held theories on graphs with Riemannian struture. 
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2.3.1. Bosonic Free Field Theory. Let X be a complex with a Riemannian struc¬ 
ture. The (massless) free held action, or Dirichlet energy, is the quadratic form 
R : C'O(X) ^ R+: 

S{f) = {df,df)cHx) 

for / e C'O(X). 

The space of helds TL C C^{X) of the discrete held theory is given by prescribing 
boundary conditions as follows: let L C dX be a nonempty subcomplex and let Ft = 
C^{X, L).The Bosonic scalar held theory is then dehned by the partition function: 

Z,{X,L)= f 
Jn 

where the measure Fijj on FL is induced by the scalar product on C^{X). The 
correlation functions are dehned in the usual way as: 

(^(Pl)---^(Pn)) = y .y [ ^(Pl) • • Vij 

Zb[A, L) 

Using standard Gaussian integral formulae give: 

log Zb{X,L) = -^logdet(AL) -h ^ (dim'H) log(27r) 

Similarly, correlation functions can be expressed in terms of matrix elements of 
For exampe, the two point function: 


2.3.2. Fermionic Free Field Theory. As before, let L C dX be a subcomplex of 
dX. Let = G°(A, L) © G^(A, L) and S~ = G^(A, L). The combinatorial Dirac 

operator ■. ^ S~ is the restriction of di + di to ., and the adjoint is the 

restriction of di + hi to S~. This operator shonld be regarded as a combinatorial 
analog of the Hodge-Dirac operator. 

The space of helds for the free fermionic held theory is Fi^ = G°(A, L)®G^(A, L) 
and Fi~ = C^{X,L). The theory is dehned by partition fnnction given by the 
Grassmanian integral over Ft = © Fi~)'. 

Zf{X, L) = exp Q(^, Dlx)^ VxV^P (2.3) 

where the measure T>x T)'ip is induced by the scalar products on and S~. Gorre- 
lation functions are once again dehned by insertions into the path integral: 

{'tl){pi)x{Pl) ■ ■ ■'>P{Pn)x{Pn)) 

= / +/ T^X 


We can explicitly calculate as follows. Ghoosing the standard basis {xi} = 
Ai\(9Ai for FL~ and {'ipj} = Ao\(9Ao U A 2 for 'H+, one can write the measure 
as: 


VxV^/j = 


n 


1 

V\Xi 
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Using the standard Grassman integral formulae, the partition function (2.4) can be 
evaluated as: 

ZKV,i) = h^^®det(£>i) (2.4) 

UVlXjl 

The fermionic held theory is closely related to the bosonic theory. First dehne 
the dual fermionic theory as: 

Zf{X, L) = exp QPli/', x)) (2.5) 

Then it is straightforward to show that: 

Zf{X,L)Zf{X,L) = {Z,{Xl,L) Z,{X,L)y^ 


2.3.3. The Dimer Model. The free fermionic theory dehned above is closely related 
to the dimer model of statistical mechanics. First we briehy recall the dehnition of 
the dimer model. On a graph G with positive edge weights w G C^{G), a perfect 
matching M of G is a subset M C Gi of edges such that each vertex v is contained 
in exactly one edge e G M. The weight of the matching M is neGM'^(®)- 
defines a Gibbs measure /i on the set A4 of all perfect matchings: the probability of 
a matching M is: 


/i(M, w) 


w{M) 

Z{G,w) 


where the partition function Z{G,w) = A gauge transformation of 

the dimer model is a function A G G°(G) that acts on the edge weights as (A • 
tc)([pi,P 2 ]) = w{\pi,p 2 ])X{pi)X{p 2 )- The gauge transformation leaves the measure 
/i invariant, but transforms the partition function as: 

Z(G,{\-w))^(ll\(p)]z(G,w) (2.6) 

VpgGo / 


The dimer model is exactly solvable by the Kasteleyn technique as follows: for 
simplicity, assume G is bipartite with black vertices {6i}j=i...n and white vertices 
{wi}j=i...n. A Kasteleyn orientation is a function w : Gi —)■ {1, —1}, such that the 
product of w on edges around any face is either —1 if the number of edges is divisible 
by four, or +1 otherwise. The Kasteleyn matrix is the n x n matrix with columns 
and rows indexed by black and white vertices respectively, and with matrix elements 

Kbiwj = w{[biWj])w{[bi,Wj]) 

when bi and Wj are adjacent, and zero otherwise. 

Kasteleyn’s theorem is that the partition function of the dimer model is given 
by Z{G,w) = det(iF). In addition, the Kasteleyn matrix encodes the statistics of 
dimer model via the local statistics theorem, which gives correlation functions of the 
dimer model in terms of the inverse Kasteleyn matrix. 

We now define a dimer model on the double of a complex X endowed with a 
local Riemannian structure. As before, let L C dX be a non-empty subgraph of the 
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boundary. Define an the edge weight function Wr G C^{D{X)) on the double as: 


WriXPa^Vr]) 


(a, a) 


Q(G) 


(d 'r)Q^i(G) 


In writing the above formula, we have assumed a is of smaller dimension than r in X. 
Let D{Xy be the graph found by attaching to each with a G dX a unit weighted 
edge and an additional univalent vertex. (See figure 2.3.3). Let Zd{X,Wr, L) be the 
dimer partition function of the dimer model defined by the graph D{Xy with weight 

Wr- 



Figure 2. A graph with boundary conditions in bold and the asso¬ 
ciated dimer model. 


Proposition 2.7. 


Za{X,Wr,L) = Zf{X,L) 


This proposition states essentially that Dc is a Kasteleyn matrix for the dimer 
model {D{Xy, Wr, L). It follows then that the correlation functions of the fermionic 
theory also give the correlation functions of the dimer model. 


Proof. First, note that the dimer model on D{Xy is equivalent to the dimer model 
found by removing from D{X) vertices Po- with a E L, and edges adjacent to them. 
In addition, applying the gauge transformation 

Jl/V^ ifaGXo,X2 

on the remaining vertices transforms the edge weights into: 


w'iiPa^Pr]) 


1 if (J G Xq, t G Xi 

|cr|/|r| if a G Xi, r G X 2 


Comparing the action of the gauge transformation (2.6) and the explicit formula for 
the fermionic partition function (2.4), it remains to show that: 

ZdiD{X)\L,w') = detiDL) 


Recall that the double D{X) has a bipartite structure. Regarding as map 
operator from white vertices in D{X) (edges in X ) to black vertices D{X) (vertices 
and faces in X we compute for t E Xp. 

(tGXo aeX2 ' ' 

where (— 1 )*^°'^) jg the function that is 0 if a is not adjacent to r, and is ±1 if the 
orientations of a and r are compatible. That is, up to a sign, the matrix element 
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is equal to the edge weight w'{\pa,pr]). It remains to check that the signs of the 
matrix elements of D satisfy the Kasteleyn orientation condition. Since every face 
in the double D{X) is quadrilateral, the Kasteleyn condition requires the number 
of minus signs around each is odd. This follows from = 0 and 5^ = 0. □ 


3. Combinatorial Laplacians in Rectangular Graphs 

In the remainder of this paper, we focus on the particular case that f/ is a rectan¬ 
gular region with a diagonal metric g = Qxx dx^+gyy dy‘^, approximated by subgraphs 
of For simplicity, we assume that both the metric and the induced Riemannian 
distance function are analytic. In addition, denote by dMi for z = 1, • • • ,4 be the 
four straight segments of dM, and let boundary conditions be specihed by L C dM, 
where L is the union of a non-empty subset of {dMi\. 

Let be the cell decomposition of M, consisting of four 0-cells placed at the 
corners of U, four 1-cells, and one 2-cell. Let {{X\ (, )c'3(xo)}iGN be a sequence of 
subdivisions by iteratively taking doubles, endowed with a local, consistent Rie¬ 
mannian structure and consistent combinatorial laplacians A^. 

Theorem 3.1. In the above setting, the determinant of the combinatorial laplacian 
A* has the following expansion: 

logdetA* = ^ f F{g) dx dy + — f ^(^f) ■ ds-F ^ log(e*)0(1) (3.1) 

C Jm JdM ^ 

where: 

F{gxx,gyy) = (^Lis + Li2 ~ Uaxxdyy) 

Bxi^dxxi dyy') — 2 ^®s (\/dxxdyy \/dxxdyy + 9yy) (3-2) 

Byi^gxx-i 9yy) 2 ^®S(\/ 9xx9yy T 9xx9yy T 9yy^ 

where Li 2 is the dilogarithm function. 


The approach to calculate the determinant is to study variations of the metric g. 
Under a variation g ^ g + Sg we have: 

4logdetA-, = Tl (3.3) 

In the next section §4, we compute an estimate of the combinatorial Green’s function 
Af, which we use in section §5 to compute an estimate of (3.3). The asymptotic 
expansion (3.1) follows by integrating (3.3), with the constant of integration hxed 
by the expansion (1.2), corresponding to the case that g = dx'^ -|- dy"^. 

We conclude this section by estimating the matrix elements of the combinatorial 
laplacians Af as z —>■ cx). 
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3.0.4. Asymptotics of the Combinatorial Laplacian. First, we fix some useful no¬ 
tation borrowed from the finite difference calculus. For clarity, we leave implicit 
superscripts indicating the level of subdivision when there is no danger of confusion, 
for example C°(X) for e for etc. For p G Xo\9Xo with coordinates 

{x,y) G M^, dehne the hnite difference operators: 

Dxf{x, y) = ^ {f{x + e,y)- f{x - e, y)) 

Dlf{x, y) = f{x + e,y)- 2f{x, y) + f{x - e, y) 


and similar operators in the ^/-directions. These operators are consistent: for every 
/ G C°°{U,dU), there are constants pf such that 


DJxf{x, y) - Ixdxfix, y) < e^pf 


2 Dlixfix, y) - Ixdlfix, y) < e^pj 


(3.4) 


Proposition 3.2. Then the laplacian has the following asymptotic asymptotic ex¬ 
pansion: 


= 


-y/ 9yy9x 


9yy jj 2 


+ 


9yy 




9yy 


(3.5) 


+o(r)Di + o(r)Di + o(€=>) A + o(€=') A 


Proof. It is straightforward to check that on the square lattice, the combinatorial 
laplacian (2.1) can be written as a sum of hnite difference operators D‘^,D‘^,Dx 
and Dy. (In fact, any operator that couples only nearest neighbor vertices and 
and annihilates constant functions can be written in terms of these hnite diherence 
operators). The proposition then follows from comparing the consistency of the 
hnite diherence operators (3.4), and consistency of the laplacian (2.2). □ 


4. Inverse of Combinatorial Laplacian 

In this section, we derive an estimate for the combinatorial Green’s function Af. 
The approach can be summarized as follows. Since the combinatorial Laplacian is a 
consistent approximation of A, the smooth Green’s function G is an approximation 
of the combinatorial Green’s function where it is smooth. However, the Green’s 
function G(p, q) is logarithmically divergent where p ^ q. We construct an approx¬ 
imate inverse matrix G® by replacing these logarithmic divergences with a lattice 
regularized logarithm function. We show that G® is a ’’parametrix,” in the sense 
that for sufficiently large i, the error 

R = I- A®G® 

tends to zero, and for sufficiently large i, the inverse A“^ can be written as 

= G\I + R + R^ + ---) (4.1) 

We begin hrst by recalling some facts about the Greens function of the smooth 
setting. 
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4.1. Smooth Green’s Function. Recall that when the Riemmanian distance func¬ 
tion d : M X M —)■ M is smooth away from the diagonal, G can be written neatly in 
terms of d. Dehne ca : M x M —)■ M as satisfying 

( 1 ) Au{pi,p2) = 0 

(2) uj{pi,p2) = log(d^(pi,p 2 )) for p2 or pi G dM 

Then the Green’s function can be written: 

^(^ 1 , 112 ) = - log {d‘^{Pi,P2)) +uj{puP2) (4.2) 

In this decomposition, both terms are divergent: the logarithmic term contains 
divergences as pi —)■ p 2 while u contains divergences as pi,p 2 dM. 

To understand the divergent behaviour as pi —)■ P 2 , one can expand the distance 
function using the Hamilton-Jacobi equation. Letting {x,y) = p 2 — Pi, we calculate 
that as (t, y) —)■ 0: 

d {pi-:P2} 9xx ^ T 9yy y 

1111 
T '^dx9xx T -|- ~dy9xx X y -\- '^dx9yy xy -\- '^dy 9 yy y 

Then for the logarithm of the distance function: 

log {d\pi,p 2 )) = log { 9 xxX^ + 9yyy‘') + log [l + 1 r* 12 ) 

\ yxx ^ “T yyy y / 

as (t, y) —)■ 0. The term on the right is bounded, so we have isolated the divergence 
into the hrst term. 

Turning to ca, recall that although for hxed pi, a;(pi, •) is smooth, as pi —)■ dM 
the function and its partial derivatives diverge. One can understand the behaviour 
of a;(pi, •) as pi —)■ dM using the method of images.(citation?). Dehne the rehected 
and the corner distance functions as: 

dR,{PuP 2 ) =inf{(J(pi,6) + d{h,p 2 ) \ b G dMi} (4.5) 

dRij{pi,P2) =mi{d{pi,bi) + d{bi,b2) + d{b2,P2) (4.6) 

I bi G dMi, 62 e dMj or bi G dMj, 62 G dMi} 

The rehected distance function is well dehned for all p 2 when d{pi, dMi) is sufficiently 
small, and similarly for the corner distance function. By dehnition, they satisfy 
dRi{pi,P2) = d{pi,p2) when pi G dMi, and dR^.{pi,p2) = dR^{pi,p2) when p 2 e dMj. 
In addition, they satisfy the Hamilton-Jacobi equation, so that A log((jR. (pi,P 2 )) = 0 
and A\og{dR..{pi,p 2 )) = 0. 

Fix S such that rehected distance function is well dehned when (Jr.(pi, •) is well 
dehned for all i and pi with d{pi,dMi). Let Bi = {p\d{p,dMi) < J} for all i, 
B = UjRj. For i j let Cij = Rj 0 Bj. Redehne u by the decomposition: 


(4.3) 


(4.4) 
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(1) For p e Cif 
G{PuP2) = 


ANANTH SRIDHAR 


- ^log(c?^(pi,P2)) + {d%{Pl^P2)) 

+ ^log {d%{Pi^P2)) - ^log (4.,(PbP2)) + a;(pi,P2) 


(2) For p e Bf. 

G{Pi,P2) =-^log (rf^(pi,P2)) +^log {d%{Pl,P2)) +Uj{pi,P2) 

(3) For p ^ B: 

G{pi,P2) = “^log {d^iPi^P^)) +^{Pi,P2) 

With this redehnition, it can be shown from a Schauder estimate that u{pi,-) is 
smooth and bounded for all pi G M, so that the divergences have been pushed into 
the reflected distance functions. 

We can compute the expansion of the reflected distance functions as before, using 
the Hamilton-Jacobi equations. Write the coordinates pj = {xj, Uj) for j = 1, 2, and 
for simplicity assume that dMi is aligned with the y-axis, and dMj is aligned with 
the x-axis. Then: 

d\^{Pl,P2) = 9xx{Pl) {xi + X 2 f + Qyyipi) {y 2 -yif + 0{x\)+0{\pi-p2\^) 
d%{pi,P2) = 9xx{pi) + X 2 f + 9yy{pi) (l/i + 92^ + O^xj + yl) + 0{\pi - PsH 
so that: 

log {d]i^{pi,P 2 )) = log { 9 ~^{xi + X 2 f + 9y^{y2 - yif) + • • • (4.7) 

log {d%^{PuP 2 )) = log { 9 x^{xi + X 2 f + 9y^{yi + y 2 f) +■■■ (4.8) 

where the remaining terms are continuous as pi,p 2 —t dM. 


4.2. Combinatorial Green’s Function. Now we construct the parametrix G* for 
the combinatorial Green’s function. It is based on lattice regularized logarithm 
function, which is a special case of a general formula given in [Ken02], arising in the 
theory of discrete complex analysis. It is dehned as follows: 


Definition 4.1. For a, b positive real numbers, dehne £^{a,b) : —)■ M for (x, y) G 

as: 


2^(a,b){x,y) = 


f log( 2 :) f z + e 


SttH 


<C Z 


hl+hl 


z + e 


hl-|: 


dz 


(4,9) 


where w = tan ^ y-^/bJaj and where G is a curve that encircles the points and 
but not the origin. 


This function satishes the following properties, the proofs of which can be found 
in [Ken02],[Buc07]: 

1 - - 2 ( 0 ,;)) ( 0 , 0 ) = 0 
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2 ' + 

3- ^{a,b){.x,y) has the asymptotic expansion for x,y ^ oo: 



ab (bx^ — (a + b)x‘^y‘^ + ay*) / 1 \ 

2471 (b x'^ + a y'^)'* + 


where 7 ^ = exp( 7 e) is the exponential of the Euler gamma 7 e. 


4. The hnite difference derivatives of il(a,b) have the asymptotic expansion: 


Dx^(a,b) {pj) y) 

^x^{a.,b) {^1 y) 

DySl/(^a,b) y) 
Dl2(^a,b){x,y) 


y/ab ax f 1 \ 

271 ax"^ + by'^ \ + 7 / 2 ) 3/2 j 

\/a6 ab y"^ — / 1 \ 

2n (ax^ + 67 / 2)2 y 2 y j 

hg +of _^_) 

271 ax2 + by^ \ (x2 + 7 / 2^372 J 

\/a 6 ab x^ — b‘^y‘^ / 1 \ 

271 (ax2 + 67/2)2 y2yj2 j 


(Here it is understood that in the hnite difference operators e = 1.) 


(4.10) 


5. At the origin, the derivative is: 


l^xii(a,fe)(0,0) =0 
L>^£(a,fe)( 0 , 0 ) = ^tan“^ 
Zl,£(,,b)( 0 , 0 ) = 0 
A’J-C(a.b)(0,0) = — tan"^ 

y 71 




4.2.1. Green’s function on X*. We now utilize the lattice logarithm dehned above to 
construct the Green’s function on X*. Roughly speaking, we replace all logarithmic 
divergences identihed in §4.1 with corresponding lattice logarithm functions. For 
some notational simplicity, we hx: 

^(Pl,P2) = - log {9xx{Pl)X^ + 9yyiP2)y‘^) 

^{P1jP 2) = ^g{pi) 

where (x, 7 /) = P 2 — Pi- We dehne also and in the obvious way 

based on the expansions (4.7). 

Then, dehne G® : Xq x Xg — )• M as: 

















14 


ANANTH SRIDHAR 


(1) For p e Cif 

G\pi,P 2) = G{pi,P2) - L{pi,p2) + £(pi,P2) 

-Lr. +£i?,(pi,P2) 

- Lr. + 2,R.{pi,P2) 

+ Lr.. - 2,R^.{pi,P2) 

(2) For p e Bp 

G\ppP 2) = G{pi,p2) - L{pi,p2) + 2^{pi,P2) 

— LR^ + £i?i(Pl,P2) 

(3) For p ^ B: 

G\pi,P2) = G{pi,P 2) - L{pi,p2) + £^{puP2) 

Here, when pi = p 2 we set Ixi{G{pi,pi) + L{pi,pi)) = 0. 

We now estimate the error i?* = A*G* — I. The series converges if in the operator 
norm |-R*|op < 1- We show the stronger condition that each column of i?* has Li 
norm of 0(e). It follows then that for sufficiently large i, \R^\op < 1- 

Proposition 4.2. Let R : Xq x Xq be defined as: 

R\Pi,P2) = A*G*(pi,P2) - dp. 

Then there is a constant C such that: 

Y, \R‘(pi,P2)\<ch 

V2&xi 


Proof. We assume that pi ^ B, and we split the sum 

\R\p,,P2)\=R\pi.Pi)R Y. (4.12) 

P2&XI P2 6AJ 

P29^Pl 

When Pi 7 ^ p 2 , we write Rfipi,p 2 ) using the fact that AG{pi,p 2 ) = 0 as follows: 

\R^\ = |A*G- A*L +A*£| 

= I A*G - A*L - e\AG - AL) + A*£ - e^AL] 

< |A*G-A*L-e^(AG-AL)| + jA^il-e^ALl (4.13) 

We estimate the two terms seperately using the aymptotic expansion of ( 3.5) of A®, 
along with the consistency properties of the hnite difference operators for the hrst 
term, and the asymptotic formulas (4.10) of A for the second term. 

We begin with the hrst term: | A®G — A®L — (AG — AL) \. Recall that the errors 
of the hrst and second order hnite diherence operators are bounded by the third and 
fourth partial derivatives of G(pi,p 2 ) and L(pi,p 2 ), which diverge near the diagonal. 
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However, we can estimate these divergent terms using the expansion of the distance 
function (4.3). Letting (x, y) = P 2 — Pi, as (x, y) —)■ 0 we have: 

G{puP2)-L{pi,p2) 




2^x9xx ^ 4“ 2^y9xx X y 2^x9yy ^9 4” 2^y9yy y 
9xxX‘^ 4“ 9yyy‘^ 


(4.14) 

(we ignore a;(pi,p 2 ) since it is smooth). From this, we can calculate explicitly that 
the third and fourth partial derivatives behave as p 2 pi as: 

. .... '"‘/in, 'n„'l=^'II -^ 

'y ^x ^y 


«“"^G(pi,P2) - d^dt-^L{p,,p2) = O 


\P2 -Plp 
1 

\P2 -PlP 


Here the constant of proportionalities in the asymptotics can be written in terms of 
the metric g and its partial derivatives, which we omit writing for brevity. Then we 
have: 

0(6^) 


\Dl{G-L)-e^dl{G-L)\ = 
\eD,{G - L) - e^d,{G - L)\ = 


\P2 -Pl\^ 
b2 -piP 


and similar expressions for the y-direction. It follows from this and equation (3.5) 
that: 

O(e^) 


|A*G'(pi,P2) - A*L(pi,P 2) - e (AG'(pi,p2) - AL(pi,p2))| = 


\P2 -PlP 


(4.15) 


We turn now to |A*£ — e^AL|. First observe that the equations (4.10) can be 
rewritten as: 


D^S.{puP 2 ) = e d^L{pi,p 2 ) + 
eD^2.{pi,P2) = e^d:cL{pi,p2) + 


0(6^) 


x^ 


0(6^) 




and similar expressions for the ?/-direction. In addition, recall that for pi ^ p 2 , we 

.iAn + sfc) i(Pi.P2)=0Thus: 


have 


Di I „2 


9 xx{pi) 9 yy{pi) ^ 

A^A-e^ALl = I (A* 


+ TidtiDl) ^ipi,P2) = 0 and 

Dl 


" ' £ - eM A 


dl 


y 


9xx{pi) 9yy{Pl) 

Using again the expression (3.5) for A®, it follows that: 


9xx{Pl ) 9yy {Pl ) 


|A®£(pi,P 2) - e^AL(pi,p2)| = 


O(e^) 


\P2 -PlP 

Putting together (4.16) and (4.15), we hnd for pi ^ p 2 '- 

O(e^) 


(4.16) 


l^*(Pl,P2)| = 


\P2 -Pl\^ 


(4.17) 
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We can estimate the sum by changing variables p = {p 2 — Pi)/^- Then: 

2 vr|p| 

b# - \ /, 

pGZ2\0 

= 0(e) 


\R\puP2)\<0{e) ^<0(e)(^^ 


P2GX5 
P2 7^Pl 


w 


-d\p\ + 1 


Lastly, we have the case pi = p 2 - We compute explicitly: 

A*£(pi,pi) = l + 0(e2) 
A*G(pi,pi) - A*(pi,pi) = 0{e^) 


(4.18) 


(4.19) 


Putting together equations (4.16) and (4.18) completes the calculation. The case 
where p E B is identical except for the additional ’’image” charges that can be 
handled in precisely the same way. 

□ 


5. Asymptotic Determinant of Laplacian 


In this section, we derive an expression (5.21) for the variation of the determinant 
of the combinatorial laplacian. It is more convenient to write variations of g in terms 
of scaling variations and . Letting A = ^jQxxQyy and uj = \/QyyjPxx-, we rewrite (3.5) 
as: 

~ + a; ^ Dy + edxUJ Dx + edyU ^ Dy^ 

+ (O(e^)D^ + 0{e^)Dy + 0{e^)Dx + 0{e^)Dy) 

Then we rewrite (3.3) as: 


Tr((Al)-^5Al) 


= E -yAUAp-Hp.p) 

p£X\dX 


- y 

X ^ 


pex\dx 


5uDl-^yDl + edx6u D, 

^ ^2 y 


dySu 
' 00^ " 


(Ai) \p,p) 

(5.1) 


+ {0{e^)Dl + 0{e^)Dl + 0{e^)Dx + 0{e^)Dy) (Al)-i(p,p) 

p£X\dX 


The hnite difference derivatives of the Green’s function G{p,p) at the diagonal can 
be estimated using the results of the last section. The key is that all relevant terms 
are given by £: 


Proposition 5.1. Asymptotics of the partial derivatives of the Green’s functions. 
For p G Ao, let rp^QM be the Euclidean distance in of p to DM. Then: 
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and similar expressions for the y-direction. When p is near the boundary, dp^oM < 
SdM, then we have a different estimate including the image charges. Let p E Bi. 
Then: 

Dl{A^)~\p,p) = Dl^(p,p) - DlS.R.{p,p) + 0{e^) 

= DAAp^p) + 0{€) 

Finally when p G Cij, we have: 

Dl{^T^{p,p) = Dl2.(p,p) - DIS.rXp,p) - DlS.R^{p,p) + Dl£.R..{p,p) + 0 ( 6 ^) 
D^{A^)-^{p,p) = -D^Ar^{p,p) - DAjip,p) + D^2 .r^.{p,p) + 0(e) 
and similar expressions for the y-direction. 


Proof. Recall that from (4.1) that: 

(A*)-^ = G* + G*R* + G*R*R* • • ■ 


We estimate 0^(A*) ^{p,p) and ^{p,p) by hrst calculating DlG’'{p,p) and 

DxG{p,p). Then we show that the higher order terms only contribute at O(e^). 

As before, we assume that p ^ B. Recall that G*(pi,p 2 ) = G{Pi,P 2 ) — ^{Pi,P 2 ) + 

2,{pi,p2). By direct computation using the expansion (4.14) for G — L, we hnd: 

DlG{p,p) - DlL{p,p) + DlA{p,p) = DlA{p,p) + 0{e‘^)D^G{p,p) - D^L{p,p) + O^A(p,p) = 0(e) 

and similar expressions for the ^/-direction. 

We turn now to the second order term: 

DlG^R\p,,p,)= Y, DlG\p,,p2)R{p2.Pi) 

Fix r to be the (euclidean) distance to the boundary, ie r = inf{|p —P 2 I | P 2 £ dM}. 

We break the sum: 

DlG\p,,p2)R{p2^Pi) = DlG\pi,pi)R\p,,pi) (5.2) 

+ DlG\p2,pi) R\puP2) (5.3) 

0 <|p 2 -Pl|<»’ 

+ DlG{P2,Pi)R\pi,P2) (5.4) 

\P2-Pl\>r 


and estimate each term seperately. 

We begin with the region |p 2 —pi\ > t’. We can estimate i?*(pi,p 2 ) as before with 
equation (4.16): 


\R\Pi,P2) \ = 


0 ( 0 ) 


\P 2 -P 1 ? 

To estimate O^G*(p 2 ,Pi), we use (4.15) and (4.10) to find: 

0 ( 0 ) 


DiG\p2,pi) = 


\P2-Pl\ 
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Thus we find: 


^lG{p2,Pi)R"{pi,P2) 

\P2-Pl\>r 


E 

\p-Pl\>r 


\P2 - Pi 




(5.5) 


Next, we turn to the region 0 < \p 2 — Pi\ < f. This region requires more work 
because there is considerable cancellation. We first find a more refined estimate for 
-R*(Pi 5 P 2 ) and show that largest terms of 0(e) in K^{pi,p 2 ) are antisymmetric about 
Pi (ie. odd as a function of {x,y) = P 2 — Pi )• We show that DlG^{pi,p 2 ), on the 
other hand is symmetric about pi. Thus, the sum (5.3) vanishes at 0(e). 

As before, we write R\pi,p 2 ) as in (4.13). Considering the term: 

|A*G(pi,P2) + A*L(pi,P2) - AG(pi,P2) - AL{pi,p2)\ (5.6) 

we can use the expansion (4.13) and expand for small {x,y): 

^ Id^gxx + ld^9yy xy"^ + ^dyPyy y^ ^ 

9xx^‘^ T 9yyy‘^ 

The first term is clearly odd in (a;, y), and the fourth partial derivatives of the higher 
order terms are only 0(l/(p2 ~ ViY) as p 2 —)• Pi- Consequently, using (3.5) again, 
we write: 

|A*G(pi,P2) + A*L(pi,P2) - AG(pi,P2) - AL(pi,P 2)| = 0(e^)/(pi,p2) + , ^ 

\P2-pir 

(5.7) 

for a function / that is anti symmetric with about pi and f{pi,P 2 ) = 

as p 2 —)■ Pi- Next, compute a similar estimate for A*£. Since both £,{pi,p 2 ) and 

L{pi,P 2 ) are even about pi it follows 

DI^{Pi,P 2) - dlL{pi,p2) 

DlQipi,P2) - dlL{pi,p2) 

are even. Since, as before we have g{pi)Dl + g{pi)Dy£, = 0 we find using (??) that 

A*X1(pi,P2) - AL(pi,P2) = e^f2{x,y) + \ (5-8) 

\P2-Pi\ 

for a function /2 that is odd and f 2 {x,y) = Putting together equations 

(5.7) and (5.8), we find: 

R\Pi,P2) = eVi(Pi,P2) + I ^ 12 (5.9) 

|P2 -Pir 

where fi{pi,P 2 ) is an odd function, bounded by . 

Next we estimate DlG^{p 2 ,Pi). Assume p 2 ^ B. First we have: 

0"G(p,,P2) - DlL(pun) = I 

bi -P 2 I 












ASYMPTOTIC DETERMINANT OF DISCRETE LAPLACE-BELTRAMI OPERATORS 19 


Turning to Dl£,{p 2 ,pi), recall that: 

Dl£.{p2,pi) = Dl£.g^p^){x/e,ye) 

We can expand out the dependence on the metric as: 

Dl2.g{P2)ixl^.yl^) = Dl£.g^p,)ip2,pi) + , I 

\Pl ~ P2\ 

Notice however that is an even function. If p 2 G B, there are additional image 
charges are bounded by ^ . In either case, we can write: 

^x<^'(Pi,P2) = e^f2(Pi,P2) + \ (5-10) 

IP 1 -P 2 I 

for a function even about pi. 

Putting the estimates (5.10) and (5.9) into (5.3), we have finally: 

^ D^G\p2,Pi)= ^ e^f2(Pi,P2)fi(Pi,P2) + = O(e^) 

' ^ ^ ^ \P 2 — Pl\ 

0 <\p 2 —pi\<r 0 <\p 2 —pi\<r 

where the hrst term vanishes due to the even-odd. 

The contribution of the higher order terms in (4.1) can be estimated using (4.16) 
and we omit the details. □ 


Returning to equation (5.1), note first that there are O(^) vertices as e —)• 0, 
so that any summand of 0{e^) in will contribute at most a constant term in the 
asymptotic expansion, which we neglect. Recall that B, C. Thus using Proposition 
(5.1), we drop irrelevant terms, and rearrange equation (5.1) as: 

TV ((Ai)-‘iAl) = 5^ [bulk] -|- [boundary] -|- [corner] (5-11) 

p£X\L peB pGC 

where the summands are: 

[bulk] = -y + j(^6uDl-^ £(pi,pi) (5.12) 

[boundary] = + edJuD^ - eS Dy^ ^Ri{Pi^P 2 ) (5-13) 

[corner] = (^^^1 - + edJujD^ - eS £^,j(pi,p 2 ) (5.14) 

we compute each of these sums in the following subsections. To compute these sums, 
we utilize the Euler-MacLaurin formula, which we write as: 


Proposition 5. 2. Let f G ([0,1]), and iV G M, then: 

Y-l , / 1 \ 

fii/N) = N f(x)dx - -(/(O) + /(!)) + O (-j 

Let S = [0,1] X [0,1] he the unit square and g G C°°{S). Then 


N-l M-1 

EE g{t/N,j/N) 

i=i j=i 


/ 9 {x,y) dx dy - 


N 


g{s)ds 0(1) 


'dS 
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5.0.2. Bulk Term. The bulk sum is: 

5X 1 / 2A ^ 5u2Xu 

—^ + T — tan (cj)---tan 

A A V TTO; Co’"' tt 

pex\L ^ 

E ^A 2 6uj / , , _i, _i.\ 

—— H-(tan(co’) — tan (u )) 

pGX\L 




Using the Euler-MacLaurin formula, we hnd that the bulk contribution to the as¬ 
ymptotic determinant is: 


f _6X 
Im ^ 

-1 r 

2e 


dM 


2 6u 

71 U 

5X 


(tan(co’) — tan ^{oj ^)) dx dy 

2 5u 


(5.15) 


71 IxJ 


(tan(co’) — tan ^(co’ ^)) ds 


5.0.3. Boundary Term. We will compute the contribution of one boundary region 
Bi. We first sum in the normal direction of the boundary c?Mj, and then sum along 
the boundary using the Euler-Maclaurin formula. We assume is aligned with the 
l/-axis. We choose coordinates (x, y) chosen so that dMi is on the y-axis with the 
interior of M in the region cc > 0. 


Observe that in the boundary summand (5.13) by symmetry the term DyS^n. 
vanishes. In addition, So we rewrite (5.13) as: 


2 Dy2R,{pi,Pi) 

We expand the terms in the cc-coordinate, normal to the boundary: 

Suj{x,y) _26uj{0,y) ^ ^ ^dJuj{0,y)X{0,y) - d^X{0,y) ^ 

\ ^ ^ A /^\0 “U 


(5.16) 


X{x,y) A(0,?/) 

dy6u{x,y) dxSu{x,0) 


mvy 


+ 


X{x,y) A(x,0) 

Recall that XlR.{pi,pi) = £,g(^x,y){2.x/e,0). We can expand the spatial dependence in 
the metric g{x,y) as: 


^xX^g{x,y){‘^x/e,0) — D^£,g(^x,o){‘2x/e) 


1 edxX{x,0) O(e^) 


+ 


(5.17) 


Stt X x^ 

This expression can be justified with using the definition of £ and the equation (4.9), 
but it is easy to guess from the asymptotic expression (4.10), which we rewrite in 
terms of A and u as: 

7^2n O(e^) 

DxX^g{x,y){^x/e,0) = ^ 

X{x,y)e‘^ 0{e^) 

DxX2g[x,y){^x/e,0) = ^ 

Using these expansions in (5.16), we hnd that all terms of order cancel so that 
we have for (5.16): 

25cu(0,?/) 2 a /o / nN , 


A(0,?/) 


'DxX2g{x,o){2x/ e, 0) 


x^ 
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We now sum normal to the boundary: 

5 






A(0,i/) 


'Dx^g{x,o){‘^x/ e, 0) + 




Su{0,y) ^ 26u{0,y) 


+ 


tan ^(cu) + 0(e) 


2^1 + a;(0,?/)2 7ra;(i/,0) 

Finally, we sum along the boundary using the Euler-Maclaurin formula: 

^ f f <5a; 2 5oj 5ijj, ^ 

- / (- , H-1-(tan ^(cn) —tan ^(o; ) di/+ 0(1) 

eJouA 2 vTTU^ tt u 7iu^ y jj ^ 

For the other boundary segments oriented along the x-axis, we have the same formula 
but with (n —We summarise this by writing: 


SB^ = 


6u 6oj 

+ — 


2vr+^ w 
6 uj Su 

OBy = - , - 

2a;\/l + uj 


And writing the boundary contribution as: 


1 /* X 

- / 6B± H —^(tan“^(a;) — tan“^(a;“^)ds + 0(1) 

e JdM 


(5.18) 

(5.19) 

(5.20) 


5.0.4. Corner Term. Using methods very similar to the previous subsection, it is 
straightforward to check that the corner terms (5.14) do not contribute to the as¬ 
ymptotic expansion. We omit the details. 


5.1. Determinant. We now compile together (5.15) and (5.20) into equation (5.11). 
The boundary term in (5.12) cancels with part of (5.13). So we have in total: 


where: 


(51ogdetA = — [ 6Fdxdy H— f 6B±ds + 0{l) 

^ JdM 


/=:2 , 
c JM 


(5.21) 


p dA 2 5cu / ,, _1, i\\ 

SF = —— H-(tan(a;) — tan ^(w )) 

X 71 OJ 


SB,, = 


Su 


2\/rT 

Su 


u^ 


Su SX 
u X 
Su SX 


SBy = , _ 

2u\/l + u'^ ^ A 

To compute the determinant, we integrate (5.21). The constant of integration is 
hxed when A = 1 and ca = 1 by the expansion (1.2). We hud: 


log det A = — / Fdx dy F 
^ J M 


B^ds + 0{1) 


'dM 
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where 

F = log(A) + [Li 2 {ioj) + Li 2 
B^ = - log(l + \/l + a;-2) - log(A) 

By = - log(l + vTTo/) - log(A) 

Substituting A = y/gxxQyy and u = \/gyyjgxx yields (3.2). 

6. Conclusions 

There are many areas where we believe this work can be improved or generalized. 
We conclude by outlining some of these areas. 

In this paper, we have focused on the square lattice on which one can approximate 
diagonal metrics. To study more general metrics, one could for example add the 
diagonal edges of the square lattice. Our methods do not generalize easily to this 
case - the difficulty is in calculating boundary contributions. In fact,the case where 
the metric is constant but nondiagonal has not been computed in the literature. One 
could also consider massive Laplacian operators of the form A + rr? or even more 
generally A + C(a;) for some potential function V G C°°{M). In this case, we expect 
the asymptotic expansion to the be same up to logarithmic terms, although we have 
not tried to prove this in this paper. For results in this direction, see [Cha06]. 

Another avenue for generalization is to study the combinatorial Laplacian on 
non-rectangular regions, for instance, rectilinear regions in We conjecture that 
the asymptotic expansion for this case is the obvious generalization of (1.3). An 
approach to this would be to derive gluing formulas for the asymptotic expansions. 
Similar gluing formulas were derived in [KenOO] by exploiting the relation to the 
spanning tree model. Alternatively, gluing formulas for determinants can be written 
in terms of the discrete Dirichlet-to-Neumann operator, which could be studied in 
the asymptotic limit using estimates for the discrete Green’s function given in this 
paper. 

In this paper, we have focused on the laplacian and have not studied the discrete 
dirac operator, and dimer model on graphs with Riemannian structure. We plan 
to study this in a future publication. Existing results already give the leading 
order contribution in the asymptotic partition function and also the average height 
function [Ken09], but we believe it is possible to strengthen these results for the 
dimer models defined in this paper. 

Finally, as we have previously mentioned, we have not attempted to investigate 
the constant term in the asymptotic expansion of the determinant. It is clear that 
the constant term depends on details of the combinatorial Laplacian not fixed by 
requiring its consistency. However, for the pure graph laplacian on the square lat¬ 
tice, corresponding to a Euclidean metric, it is known that the constant term in 
the expansion is equal to the C-regulai'ized determinant of the analytic laplacian. 
We question if there is a natural condition on combinatorial Laplacians for which 
the constant term in the asymptotic expansion of the determinant is equal to the 
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(■-regularized determinant. Identifying such a condition would be a step towards 
defining a free Gaussian field theory as the regularized limit of discrete Gaussian 
held theories. 
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